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The scaling of the slow structural relaxation with the fast caged dynamics is evidenced in 
the molten salt Cao.4K o.6(NO 3) 1.4 (CKN) over about thirteen decades of the structural 
relaxation time. Glycerol scaling was analyzed in detail. In glycerol, the short-time mean- 
square displacement («^), a measure of the caged dynamics, is contributed by free- volume. It 
is seen that, in order to evidence the scaling, the observation time of the fast dynamics must 
be shorter than the time scales of the relaxation processes. Systems with both negligible (like 
CKN, glycerol and network glassformers) and high (like van der Waals liquids and polymers) 
pressure-energy correlations exhibit the scaling between the slow relaxation and the fast caged 
dynamics. According to the available experiments, an isomorph-invariant expression of the 
master curve of the scaled data is not distinguishable from a simpler not-invariant expression. 
Instead, the latter agrees better with the simulations on a wide class of model polymers. 
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1. Introduction 

Understanding the extraordinary viscous slow-down that accompanies glass forma- 
tion is a major scientific challenge On approaching the glass transition (GT), 
trapping effects are more and more prominent. The average escape time from the 
cage of the first neighbors, i.e. the structural relaxation time Tq,, increases from a 
few picoseconds up to thousands of seconds. The rattling motion inside the cage 
occurs on picosecond time scales with amplitude (u^)^/^. This quantity is related 
to the Debye- Waller factor which, assuming harmonicity of thermal motion, takes 
the form exp (^—q'^{u'^)/3) where q is the absolute value of the scattering vector. 
At first sight, due to the extreme time-scale separation between the rattling mo- 
tion (~ 10~^^ s) and the relaxation (tq, ~ lO^s at GT), one expects the complete 
independence of the two motions. Nonetheless, several authors investigated their 
correlations, emphasizing in particular the link with the bulk elastic properties 
(for a review see ref. In this research field, the universal scaling between the 
structural relaxation time, or the shear viscosity, and (u^) was reported for several 
numerical models, including linear polymers, mixtures, prototypical glassformers 
like Si02 and o-terphenyl (OTP), and icosahedral glassformer [sl-uj. The resulting 
master curve fits with the available experimental data from supercooled liquids, 
polymers and metallic glasses over about eighteen decades of relaxation times (or 
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viscosity in non-polymeric systems) and a very wide range of fragilities. 

The present paper shows that the scaling between the structural relaxation time 
and (n^) also holds for the molten salt Cao. 4X0.6(^03) 1.4 (CKN) over about thir- 
teen decades of Ta- Glycerol scaling was analyzed in detail. The scaling of glycerol 
and, over a less wide range of relaxation times and viscosities, network glassform- 
ers like Si02, Ge02 and B2O3 was already reported 0-0] • Here, for the first time, 
the scaling of these systems and CKN will be discussed from the viewpoint of the 
pressure-energy correlations . In addition, it is shown that the role of the free 

volume in glycerol is not negligible and that the scaling works only if the obser- 
vation time of the fast dynamics is shorter than the time scales of the relaxation 
processes 0]. 

The paper first summarizes the results about the scaling between the structural 
relaxation and (u^). Then, the results on glycerol and CKN are discussed in detail. 
Finally, the conclusions are presented. 



2. Scaling between viscous flow, structural relaxation and (u^) 

Fig(T] shows the scaling between the structural relaxation time and the viscosity 
7] with respect to the reduced mean-square displacement (MSD), {u'^)/{ug), with 
(tig) = {v?{Tg)). The experimental data include supercooled liquids, polymers, 
metallic glasses over about eighteen decades of relaxation times and a very wide 
range of fragilities. The plot now includes the ionic liquid CKN as well. All the 
data in FigH] are at ambient pressure. The master curve of FigH] (black line) is 
expressed analytically by: 

with X equal to the reduced quantities Tq/to or rj/rjQ. The best-fit values ( 
a = —0.424(1), P = 1.62(6) and 7 = 12.3(1) ) were drawn by Molecular-Dynamics 
simulations (MD) on model polymeric systems and mixtures, being confirmed 
by comparison withprototypical glassformers like Si02 and o-terphenyl, and icosa- 
hedral glassformer [6|. For the systems in Fig(T] Table [T] lists the MSD at the glass 
transition temperature, i.e. the temperature where = 10^ s or 77 = 10^^ Pa • s, 
and the conversion factors tq and 770 between the actual time and viscosity units 
and the corresponding MD units, respectively. Note that the conversion factors are 
the only adjustable parameters of the overall scaling procedure. Table [1] shows also 
the (approximate) observation time of the fast dynamics by the experiments At. 
The structural relaxation times of the systems in Fig. [1] are comparable or longer 
than At. The At dependence of the scaling of glycerol is discussed in Sec. El 



3. Scaling and the pressure-energy correlations 

Figl2] shows that, when the CKN structural relaxation time and the glycerol shear 
viscosity are plotted vs the reduced MSD, they collapse over about thirteen decades 
in a single master curve well described by Eq{T] within the error bars. 

There is competition between van der Waals and Coulombic terms in the inter- 
acting potential of CKN and glycerol. This feature deserves consideration from the 
viewpoint of the pressure-energy correlations Such correlations are expected 

to be weak for CKN and glycerol, thus leading to the absence of isomorphic states 
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in these two systems [41[. Weak correlations are also expected for network glass- 
formers like Si02, Ge02 and B2O3 which exhibit scaling (see FiglH not included 
in Fig 12] for clarity reasons) [SkZI]. 

Strongly correlating systems, e.g. van der Waals liquids, have isomorphic states. 
Rigorously, any general theory of the liquid state must end up in relations expressed 
only in terms of isomorph invariants (constant quantities when evaluated over a set 
of isomorphic states) to deal with the strongly correlating systems (a criterion also 
known as the "isomorph filter" [l3])- From this respect, the master curve given by 
EqU] with constant a, /3,^ parameters does not pass the "isomorph filter" and then 
should be unable to encompass strongly correlating systems since: 

where C = C/Co with ^0 equal to p~^^^, p~^/^ ^ym/ksT, p^/^^JmkBT for length, time 
and viscosity, respectively and X = t^^/tq, fj/ffo [4l[. A proper isomorph- invariant 
modification of Eq[T]is: 

,„,X = C.+fti|| + C,(^i^j (3) 

where Ci, i = 1,2,3 are constants. The above equation is more complex in nature 
than Eq.[T] since it needs the density as additional input parameter. 

The fact that Eq. [3]is constrained by the basic properties of strongly correlating 
systems, which are not accounted for by Eq. (H poses the question of the flexibility 
of these master curves to deal with systems with either strong (like OTP, TNB or 
polymers) or weak (like CKN, glycerol and the network glassformers Si02, Ge02 
and B2O3) pressure-energy correlations. May Eqs. [T] and [3] be discriminated? To 
date, little help is provided by the experiments. The experimental data plotted 
in Fig{T] are taken at ambient pressure by sweeping the temperature in a limited 
range, thus resulting in small density changes. Owing to the much larger changes 
of the structural relaxation time, the viscosity and MSD, the equalities logX ~ 
logX, (u^) / (v?) ~ {Ug)/{u^) hold to a good approximation and Eqs. [H [3] are 
hardly distinguishable. To our knowledge, joint data concerning Tq, (or viscosity for 
non polymeric systems) and MSD from high-pressure experiments spanning larger 
density changes are not available yet. If no clear-cut conclusion on the experimental 
side may be still drawn, numerical studies offer more insight. In particular, Eqs.[T] 
and [3] were investigated in the temperature/density phase space of a class of model 
polymers with different chain lengths and generalized Lennard-Jones potentials 
having constant position and depth of their minimum 0. Note that pressure- 
energy correlations change with the potential and the physical state [44]. The 
different systems were compared with no MSD or Tq, rescaling and Eq(T]was used 
as master curve by taking (Ug) as a constant independent of the model polymer. 
These positions: i) rely on the finding that the raw curves vs (u^) for different 
systems superimpose very well, ii) fit in with the need to limit the number of 
adjustable parameters. To keep on an equal footing the comparison between the 
two master curves, an analogous position was adopted for Eq. [Sj i.e. {Ug) was set to 
a constant independent of the model polymer (a step still preserving the isomorph 
invariance of Eql3]). On this basis, Eq. [T]was found to be more consistent with the 
numerical data (see Fig. 11 of ref. [a]). To summarize, if the experiments do not 
discriminate between the two versions of the master curve, the isomorph-invariant 
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version Eq. [3l which is more complex in nature than Eq. [T] as noted above, agrees 
less with the simulations of a wide class of polymeric systems. 



4. Free volume efTects in Glycerol 



The closeness of (u ) with free- volume concepts was noted in experimen ts |43l .l46- 
4^ . theories 0, [4§, 5C] and simulations 0,0, 51] with some debate 52. l53l]. The 
basic question is about the existence of some critical displacement above which 
structural changes take place and it was stated that the rearrangement only takes 



place when there is a local, temporary density decreas e [4 ]. We refer also to previous 



experimental [SJ, l55|] , theoretical and simulation [5614591] work. 

The ratio between the volume that is accessible to the monomer center-of-mass 
and the monomer volume is Vq^^^^ ~ From our 

simulation data (first line of Table [1]) one finds: 

4^^^) ~ 0.017 (4) 

Flory and coworkers proposed that the glass transition takes place under iso-free 
volume conditions with the universal value vq ~ 0.025 [g^]. This supports the 
conclusion that (u^) and the free volume are related, as above noted. 

To better elucidate the matter, we consider results from the Positron Annihi- 
lation Lifetime Spectroscopy (PALS) [H, [H, [i^. PALS is a structural technique 
to parametrize the free volume. In fact, the orthopositronium bound state of a 
positron has a strong tendency to localize in holes of low electron density and de- 
cays with lifetime being related to the average cavity size. In particular, we test 
the ansatz : 

{n^)=Cr3 (5) 

where C is a constant. The above equation captures the essentials of the relation 
between rs and («^), both increasing with the free volume. Furthermore, it provides 
a simple procedure to compare PALS and (u^) data with no adjustable parameters. 
In fact, Eqj5] yields the equality {v?g) / {v?) = T-sg/r-s, with r^g = T^iTg)^ i.e. the 
reduced MSD is equal to the reduced PALS lifetime t^. In Fig. [3] the viscosity data 
of glycerol were plotted by replacing the reduced MSD with the reduced rs and 
comparing the results to the master curve expressed by Eq.[TJ The conversion factor 
between the actual and the MD viscosity units was adjusted to logr^g = logryo — 1. 
Other parameters like in Table [H Fig. [3] shows that the scaling of the long-time 
dynamics by either the reduced rs or the reduced MSD are quite close to each 
other over about ten decades in relaxation times, thus leading to the conclusion 
that {u^) and free- volume are correlated. Deviations are observed for glycerol when 
cage restructuring, and then the free- volume fluctuations, occurs on time scales 
shorter than the PALS observation time, i.e. Tq, < rs ~ 2 ns (see also ref.fi^]). In 
this regime the cage dynamics still takes place (it is detectable for r^ > 1 — 10 ps) 
and {v?) is still well-defined (1,B] but the observation time rs is too long and almost 
temperature- independent [43|] . 

It must be noted that PALS data are usually interpreted by a model assuming 
the spherical shape of the hole where the positron decays. The model relates rs to 
the hole radius R [13] . In addition to EqlSj we also tested the relation (u^) oc {R 
drawn by the fit of rs data) and found that the correlation plot between log r/ and 
R^ agrees with the master curve expressed by Eq. [1] with 5^1. However, being 
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a model-dependent conclusion involving adjustable parameters, we think that the 
analysis of the free- volume role in terms of Eq. [5] is more robust. 

5. Time scales of (it'^) in Glycerol 

We now proceed and better clarify that the correct investigation of the cage rattling 
is ensured if the observation time of the fast dynamics is shorter than Tq. Fig. [J] 
shows the correlation between (u^) and the viscosity data of glycerol, (u^) was 
taken as the mean square displacement within the time At, the latter depending 
on the energy resolution of the neutron scattering (NS) experiment. It is seen that 
the scaling between the long-time dynamics and {v?) At=OAns is rather good even 
for Ta/ /S.t ~ 1, whereas using {u'^)At=5ns leads to deviations from the universal 
master curve if r^/At < 60. This suggests that, if Ta ^ 300 ns, {u^)At=5ns is 
contributed by components which are not related to the cage vibrational dynamics. 
One anticipates that these components have relaxational character. Interestingly, 
T3 from PALS tracks the universal master curve of the structural relaxation down 
to Tq ~ 3 ns where Ta/ AtpALS ~ Tq/ts ~ 1 (Figj3]). The finding seems to indicate 
that in the range 3 ns < < 300 ns the free-volume is less affected by the 
relaxation than MSD on the nanosecond time scale. The above discussion supports 
the conclusion that cage-rattling MSD must be measured on time scales At shorter 
than the time scales due to both the structural and the possible local relaxations. 
Shortening At results in a remarkable extension of the region where the scaling 
between the structural relaxation and the cage dynamics of glycerol is observed 
(Fig. |4]). However, this effect is less apparent in OTP 

6. Conclusions 

It is shown that one molten salt (CKN) and one hydrogen-bonded liquid (glycerol) 
exhibit the scaling between the mean-square displacement due to cage rattling, 
(u^), and the structural relaxation time or viscosity already observed in super- 
cooled liquids, polymers and metallic glasses. Systems with both negligible (like 
CKN, glycerol, Si02, Ge02 and B2O3) and high pressure-energy correlations (like 
OTP, TNB or polymers) exhibit the above scaling. According to the available ex- 
periments, an isomorph-invariant expression of the master curve of the scaled data 
is not distinguishable from a simpler not-invariant expression. Instead, the latter 
agrees better with the simulations of a wide class of model polymers. The cage 
rattling of glycerol is contributed by free-volume over about ten decades in relax- 
ation times. Both PALS and neutron scattering experiments show that, in order 
to evidence the scaling, the observation time of the fast dynamics must be shorter 
than the time scales of the relaxation processes. 
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TABLES 

Table 1.: Details about the systems in Fig(T] (arranged in order of in- 
creasing fragility) and the MD simulations used to derive Eq(TJ The 
experimental structural relaxation time Tq is drawn by dielectric 
spectroscopy apart from B2O3 and CKN. tq and r/o ^-re the conver- 
sion factors between the actual time and viscosity units and the cor- 
responding MD units, respectively. Note that, apart from B2O3 and 
CKN, log To and log 770 cover the narrow ranges — 1 > logro > —2 
and —11 > logr/o > —12. MSD is drawn by Incoherent Neutron 
Scattering (INS) or Mossbauer Spectroscopy (MS). The table lists 
the approximate observation time of the fast caged dynamics by the 
experiment At, the MSD at the glass transition (u^) (in A ) or, 
equivalently, the Lamb- Mossbauer factor — In/^. 



System 










MSD 


















quantity 


logTo.logrjo ref. 


technique (At) 


(«g),-ln/g ref. 


MD 









MD 


0.01667 




SiOa 


V 




If 


INS (40ps) 


0.081 


if 


Ge02 


r] 


-1 


[131 


INS (0.4ns) 


0.191 


[141 


B2O3 


V 


+2.2 


[151 


INS (40ps) 


0.065 


[161 


B2O3 


rj 


-8.4 


[17] 


INS (40ps) 


0.065 


[16] 


V4t alloy 


V 


-1 


[181 


MS (O.lAis) 


0.885 


[191 


Glycerol 


V 


-1 


\20\ 


INS (0.4ns) 


0.022 


pTl 


1,4 PI 




-12 


p2l 


INS (4ns) 


0.427 


p31 


TNB 


V 


-2 


\Ml 


INS (0.4ns) 


0.315 


p51 


Fo+DBP 




-11 


[261 


MS (O.l^is) 


3.15 


1261 


Fc+DBP 


V 


-2 


[271 


MS (O.l^is) 


3.05 


[261 


OTP 


Ta 


-11 


125] 


INS (0.4ns) 


0.215 


psi 


OTP 


V 


-1 


[281 


INS (0.4ns) 


0.232 


[281 


Selenium 


V 


-1.66 


[29] 


INS (20ps) 


0.155 


[29] 


CKN 


tJ 


-13.5 


[30] 


INS (lOps) 


0.132 


[31. 321 


1,4 PBD 


Ta 


-11 


[33] 


INS (4ns) 


0.102 




a-PP 


Ta 


-11.5 


[361 


INS (0.2ns) 


0.13 


[371 


PMMA 


Ta 


-11.5 


[381 


INS (5ns) 


1.1 


1391 


PVC 


Ta 


-11 


[40] 


INS (5ns) 


0.51 


[39] 



Zr46.8Ti8.2Cu7.5NiioBe27.5 

Data aggregated from different techniques. 
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Figure captions 



Figure 1.: Reduced relaxation time and viscosity vs reduced MSD factor {{v?g) = 
{v?iTg))). Tlie numbers in parenthesis denote the fragility m. The black curve is 
Eq{TJ The colored curves bound the accuracy of EqU] d, B]. Refer to Table [D for 
further details about the experiments. 



Figure 2.: Reduced relaxation time (CKN) and viscosity (glycerol) vs the reduced 
MSD. The other curves have the same meaning of Fig. [TJ The numbers in paren- 
thesis denote the fragility m. Other parameters are listed in Tabled) 



Figure 3.: Test of Eq. p for glycerol by using two different data sets for , i.e. 
ref. [4^ (Glyl) and ref.jS] (Gly2). The plots of the reduced glycerol viscosity vs. 
the reduced time ts/ts^ ( Tsg = T3{Tg) ) are compared with the master curve EqU] 
(black curve). The other curves have the same meaning of Fig. [H tq and tjq are 
given in Tabled! logr/Q = logr/o — 1. Deviations are observed when the free- volume 
fluctuations become fast with respect to the PALS timescale, i.e. < ~ 2 ns. 



Figure 4.: Reduced glycerol viscosity vs. the reduced MSD taken at two different 
timescales At = h/AE depending on the energy resolution of the neutron scatter- 
ing experiment AE ( IN13 and IN16 data from refs.[2ll] and 45|], respectively). All 
the curves have the same meaning of Fig. dl tq and rjQ are given in Table dl The 
arrows point to states with the indicated values. 
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